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EQUATION WITH MULTIPLICATIVE NOISE 

By Anne de Bouard and Arnaud Debussche 



' We study the influence of a multiplicative Gaussian noise, white 

^N^j . in time and correlated in space, on the blow-up phenomenon in the 

supercritical nonlinear Schrodinger equation. We prove that any suf- 
ficiently regular and localized deterministic initial data gives rise to 
Ph a solution which blows up in arbitrarily small time with a positive 

Qh ' probability. 

■s 

d ■ 1. Introduction. The understanding of the influence of a noise on the 

propagation of waves is a very important problem. Although the propagation 
is often described by deterministic models, in many circumstances random- 
ness should be taken into account. It can change drastically the qualitative 
J> \ behavior and result in new properties. 

■ This is a very vast subject. Propagation can be described by many differ- 

(y-j , ent models and randomness can take several different forms. In this article, 

' we consider the nonlinear Schrodinger equation as the deterministic model. 

It describes the propagation of waves in media with both nonlinear and 
dispersive responses. It is used in many areas of physics, for example, hy- 
drodynamics, plasma physics, nonlinear optics, molecular biology, and so 
on. 

It is well known that this equation has localized solutions called solitary 
waves or, sometimes, solitons. When the nonlinearity is not too strong — or 
subcritical, these are particularly robust and propagate without changing 
form. However, for stronger nonlinearity, these are unstable and the insta- 
bility results in the collapse — or blow-up — of the localized wave. 

Here, we are particularly interested in the influence of a noise acting as a 
potential on this behavior. Such noise has been considered in [14]; there the 
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paths of the noise are smooth functions and the nonlinearity is subcritical. 
In the case of a white noise, which is considered here, this type of model 
has been introduced in the context of crystals (see [1, 2] and also [18, 21] 
for other models). It is expected that such a noise has a strong influence on 
the solutions which blow-up. It may delay or even prevent the formation of 
a singularity. In [12], some numerical simulations tend to show that this is 
the case for a very irregular noise: for a space-time white noise. However, 
in the supercritical case and for a noise which is correlated in space but 
nondegenerate, it has been observed that, on the contrary, any solution 
seems to blow-up in a finite time. Recall that in the deterministic case, only 
a restricted class of solutions blow-up. Our aim is to prove rigorously such 
a behavior. 

It is mathematically very difficult to consider space-time white noises; this 
is due to the lack of smoothing effect in the Schrodinger equation. Thus, we 
restrict our attention to the study of correlated noises. 

The case of an additive noise has been considered in [8, 9] and it has 
been proved that for any initial data, blow-up occurs in the sense that, for 
arbitrary t > 0, the probability that the solution blows up before the time 
t is strictly positive (see Remark 4.4). Thus, the noise strongly influences 
this blow-up phenomenon. This result is in perfect agreement with the nu- 
merical simulations. The argument is based on three ingredients. First, we 
generalize the deterministic argument to prove that blow-up occurs for some 
initial data. This is based on a stochastic version of the variance identity 
(see [20, 22]). Then, we use the fact that the nonlinear Schrodinger equation 
is controlable by a forcing term. Thus, any initial data can be transformed 
into a state which yields a singular solution. Finally, since the noise is non- 
degenerate and the solution depends continuously on the path of the noise, 
we can argue that, with positive probability, the noise will be close to the 
control so that blow-up will happen afterward. 

In the multiplicative case, we can again generalize the variance identity 
and prove that a restricted class of initial data evolves into a singular solu- 
tion. An initial data with sufficiently negative energy is in this class. This 
is done in Sections 3 and 4, after we have set some notation and recalled 
some preliminary results in Section 2. However, it is not known if the nonlin- 
ear Schrodinger equation is controlable by a potential and it is well known 
that the solution does not depend continuously on the noise. Thus, it is 
not straightforward to generalize the argument used in the additive case. 
In Section 5 we address this problem and show that, using more subtle ar- 
guments, we can prove the same result as in the additive case, when the 
nonlinearity is supercritical. First, for any initial data, we construct a de- 
terministic potential such that, when it is used as a control in the nonlinear 
Schrodinger equation, the energy of the solution becomes arbitrarily nega- 
tive. The problem is then to replace the continuity argument. In fact, we 



STOCHASTIC NONLINEAR SCHRODINGER EQUATION 



3 



prove that the controlled solution is in the support of the diffusion associ- 
ated to the stochastic nonlinear equation, which is exactly what we need to 
prove our result. The techniques used to obtain this last step are inspired 
from [19] and we refer the reader to this article for a nice and simple presen- 
tation of support theorems, as well as for references on this subject. Many 
ingredients used in [19] are similar to the ones used to prove the convergence 
of a numerical scheme. Thus, we have been led to use the ideas developed 
in [11]. 

2. Notation. We consider a complete probability space (Q,J-,F), en- 
dowed with a filtration (Ft)t>0i an d a sequence (Pk)k&N of independent 
real- valued Brownian motions on 1R + associated to the filtration (J-t)t>o- 
We study NLS equations with multiplicative noise. The noise is real valued. 
To define it rigorously, we use the space L 2 (M d ;M) of real- valued square m- 
tegrable functions on M. d and a Hilbertian basis (ek)k& of this space. Then, 
given a bounded linear operator (p on this space, we define the process 

oo 

W(t,x,u) = Y J Pk(t,u)<f>e k (x), t>0, x£R d , weO. 

k=0 

It is a Wiener process on L 2 (R d ;M), with covariance operator (fxfi* . In all 
that follows, <p4>* will be assumed to be a finite trace operator in L 2 (]R rf ;]R); 
this easily implies, for a fixed t, the convergence of the series above in 
L 2 (^;L 2 (M d ;M)), and almost surely in L 2 (R d ;R). The time derivative of 
W models a noise which is delta correlated in time. If <fi is defined through 
a kernel /C, which means that for any square integrable function u, 

cf)u(x) = I JC(x,y)u(y)dy, 
Jw. d 

then the correlation function of the noise is formally given by 
fdW dW \ 

with 

c(x, y)= lC(x, z))C(y, z) dz. 

Let a > 0, we consider the following stochastic NLS equation, introduced in 
[1] and [2], 

(2.1) idu - (Au + \u\ 2a u) dt = u o dW, 

where o stands for a Stratonovich product in the right-hand side of (2.1). 
We will use the equivalent Ito form of this equation. Define, for x € M. d , the 
function 

oo 

F <t>( x ) = J2(4>ek(x)) 2 , 

k=0 
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which does not depend on the basis (ek)keN] this equivalent ltd equation 
may be written as 

(2.2) idu - (At* + \u\ 2a u) dt = udW - ^uF^ dt 

(see [7] for details). 

A local existence result for (2.2) has been proved in [10]. Before recalling 
it, we need to set the functional framework. 

For p > 1, L p (R d ) is the classical Lebesgue space (of complex valued func- 
tions), and the inner product in the real Hilbert space L 2 (R d ) is denoted 
by (•, •), that is, 

(u,v)=~Re / u(x)v(x)dx 



for u,v 6 L 2 (R d ). The norm in L p (R d ) is denoted by | • \lp and, for p = 2, 
by | • |- 



We define the usual space H S (R ) of tempered distributions v £ S'( 
whose Fourier transform v satisfies (1 + \£\ 2 ) s/2 v € L 2 (R d ); H s (R d ,R) is 
the subspace of H s (R d ) consisting of real- valued functions. The norm in 
H s (R d ) is denoted by \-\h s - For p e N, W 1,p (R d ) is the space of functions in 
L p (R d ) whose first order derivatives are in L p (R d ) and its norm is denoted 
by I • Iw 1 ^- 

If / is an interval of R, E is a Banach space, and 1 < r < +oo, then 
L r (I;E) is the space of strongly Lebesgue measurable functions v from I 
into E such that the function 1 1— > \v(t)\E is in L r (I). We define similarly the 
spaces L r (Q;E). 

Given two separable Hilbert spaces H and H, we denote by C2(H, H) the 
space of Hilbert-Schmidt operators $ from H into H, endowed with the 
norm 

|<3>| 2 ~ = tr$*$ = ^ l$ei-l~ 
fceN 

where (efc)fceN is any orthonormal basis of H. When H = L 2 (R. d ;M.) and 
H = H s (R d ;R), then C 2 {H,H) is simply denoted by £% s . Given a Banach 
space -B, we will also consider bounded linear operators from L 2 (R d ) into 
B, and in order to replace the notion of Hilbert-Schmidt operators, we 
use in this case the notion of 7-radonifying operators (see, e.g., [3, 4]). We 
denote by R(L 2 ,B) the space of 7-radonifying operators from L (R ;R) 
into B. We recall that if B is a Hilbert space, then R(L 2 ,B) coincides with 
£ 2 (L 2 {R d ;R),B). 

It is classical that if ip(-) is a predictable random process defined on [0, T], 
with values in the space of continuous operators from L 2 (IR d ,IR) into an 
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Hilbert space K, and is such that 

i-T 







( s ) ( / ) \ 2 c 2 (L 2 (R d ,R),K) ds <oo a.s., 



then the stochastic integral 



e^JO 



(el 

is a well-defined i^-valued random variable. 

We now recall the existence result of [10]. For technical reasons, we restrict 
our attention to the case d < 3. However, if the assumptions are suitably 
modified, most of the results can be extended to larger spatial dimensions. 

We recall that a pair (r, p) of positive numbers is called an admissible pair 
if r>2 and | = d(i-|). 

Theorem 2.1. Assume that < a if d=l or 2 and < a < 2 if d = 3. 
Let 4> G £2' H R(L 2 ,W 1,a (M. d )) with a > 2d. Then there exists an admissible 
pair (r,p) such that for any To -measurable uo with values in H l (M. d ), there 
exist a stopping time t*{uq,lv) and a unique solution to (2.2) starting from 
u which is almost surely in C([0,r]]H 1 (R d ))nL r {0,T]W 1 ^(R i )) for any 
t < t*(uo). Moreover, we have almost surely, 

t*(uo,lo) = +oo or limsup \u(t)\ H i(^d\ = +oo. 

tyr*(uo,oj) 

In this article we want to study blow-up phenomena. As in the determin- 
istic case, we need to work in the context of H 1 solutions. This is why we 
use the result of [10]. In [7], L 2 solutions are considered, requiring weaker 
assumptions on (ft; however, there, a stronger restriction is imposed on a and 
it is not expected that blow-up occurs in that case. Indeed, it is shown in 
[7] that the solutions are global in time. In fact, we will require that a does 
not fulfill the assumption of [7]. We will work in the so-called critical and 
supercritical cases: 

2 

The critical case corresponds to a = I and most of our results will not be 
true in that case. Note that if a satisfies the assumptions of Theorem 2.1, by 
the Sobolev embedding theorem, we know that i? 1 (]R ci ) C L 2a+2 (M. d ) with a 
continuous embedding. 

In fact, we need stronger assumptions on the initial data. As in most of 
the deterministic blow-up results, we require some spatial localization. For 
rj > 0, we introduce the spaces 

S 7 ' = {v G iF»(R d ) : [x|"« G L 2 (R d )} 
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endowed with the norm 



I I 2 I |2 ,|| m |2 



When r\ = 1, we set S 1 = S. 

In all of the article, c or c(-, ...,•) is a constant which may change from 
one line to another and depends only on its arguments. 

3. The stochastic variance identity. In this section we derive an identity 
on the evolution of the variance, or virial, of the solutions of the nonlin- 
ear Schrodinger equation (2.2). This is a generalization of the well-known 
corresponding formula in the deterministic case. 

The variance (this quantity should not be confused with the probabilistic 
variance) is defined by 

V(v)=[ \x\ 2 \v(x)\ 2 dx, ueS. 

Its evolution is described in terms of the mass, 

M(v)= f \v{x)\ 2 dx = M 2 , veL 2 (R d ), 

jR d 

the energy, 

H{v) = \f \Vv(x)\ 2 dx - f \v(x)\ 2 ° +2 dx 
2 JwA 2a + 2 

and the momentum, 

G(v) =lm v(x) x ■ Vv(x) dx, v 6 S. 

The first result describes the evolution of M and H and is proved in [10]. 
We introduce the function f \ depending on the covariance operator and 
such that, for any orthonormal basis (eg)i^ of L 2 (M d ;]R), 

fl(x)=J2\V<pe e (x)\ 2 , x£R d . 

Note that this function does not depend on the chosen basis (e^gpj. Indeed, 
since </> is assumed to be Hilbert-Schmidt from L 2 (R d ) into H l (R d ), it is 
associated to a kernel K € L 2 (R d x R d ) such that V ' X K € L 2 (R d x R d ); it can 
be seen that f^(x) = \V x fC(x, -)| 2 2{ 
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Proposition 3.1. Let uq, a and (j) be as in Theorem 2.1. For any stop- 
ping time t such that r < t*(uq) a.s., we have 

(3.1) 

and 



M(u(t)) = M(u ) a.s. 

H(u(t)) = H(uq) — ImY^ / / u(s, x)Vu(s, x) ■ V(^e^)(x) dxd(3i{s) 

fztJo Jm. d 



(3.2) 



u(s , x)\ 2 f^(x) dxds a.s., 



o JM. d 



where u is the solution of (2.2) given by Theorem 2.1 with u(0) = uq. 

The different terms in (3.2) make sense. To see this we first observe that, 
since </> € C 2 ' j the function is an L 1 function and since 4> € R(L 2 ,W 1,a (M. d )), 
we also have f\ € L a / 2 (R d ) (see [3, 4, 7]). Thus, thanks to Holder's inequal- 
ity, /I € L 9 (]R rf ) for any q 6 [l,a/2]. We then introduce the stopping time: 

T fe = inf{se [0,t], |«(s)|^i >&}• 
For any fc£N, we have 

2\ 



E 



E 



E 



w(s, x)Vu(s, x) ■ V((/>e^)(x) dxdfi^s) 



VfeN 



TATfc 



u(s, x)V!i(s, x) • V((^e£)(x) tix 



< E 



rAr fe 



|u(s, x)\ \Vu(s, x)\(f^,(x)) 1 ^ 2 dx j ds 



thanks to Minkowski's inequality in the last step. Using Sobolev's embed- 
ding, it is now easy to see that this is a finite quantity. Since, for almost 
every we have r < for some k, it follows that the stochastic term 

in (3.2) is a well-defined stochastic integral. 
We now investigate the evolution of V and G. 

Proposition 3.2. Let the assumptions of Proposition 3.1 hold and as- 
sume, moreover, that uq £ £ a.s., then for any stopping time r such that 
t<t*(uq) a.s., we have 

G{u(r)) = G(u ) + A T H{u(s))ds+ 2 



(3.3) 



o- + l Jo 



Hs)\ 2 °+l2 ds 



+ E 

fceN 



JM. d 



\u(s,x)\ 2 x ■ V{4>ek){x) dxdPk(s) 
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and 

(3.4) V{u{t)) = V(u ) + 4 r G(u(s)) ds. 

Jo 

Remark 3.3. The proof of this result is postponed to Section 6.1. Note 
that, formally, it is the result of the same computation as in the deterministic 
case and that it is easier to perform the computation with the Stratonovich 
equation since in that case the standard calculus can be used. However, in 
the proof we use the Ito form which is better suited for rigorous justifications. 

A similar result was obtained in the case of an additive noise in [8] . The 
formula was more complicated there. In the case of a multiplicative noise, 
more terms cancel in the computation. This strongly uses the fact that the 
noise is real valued and reflects the fact that a Stratonovich multiplicative 
noise corresponds to a noisy potential and is a more realistic representation 
of physical phenomena in the present context. 

Remark 3.4. We will actually prove that, for each t > and k € N, there 
exists a constant c(V(uo),t,k) such that V(u(t)) < c(V(uq), t, k) provided 
t<T*(u ) and \u(s)\ H i < k for s £ [0,t] [see (6.2)]. 

Using this and similar arguments as above, it can be seen that the stochas- 
tic term in (3.3) is well defined. 

Corollary 3.5. Under the same assumptions as in Proposition 3.2, 
we have 

V(u ) + 4G(u )t + 8H(u )t 2 

+ 4 ^T f f Hsi)\%tld Sl ds 
a + 1 jo Jo 

pT rs ps\ p 

+ 8/ / / / \u(s2,x)\ 2 fl(x)dxds2dsids 
Jo Jo Jo JR d 

rr rs r 

+ 4 Y1 / / / \u(si,x)\ 2 x- V{4>e k ){x)dxdf5 k {s 1 )ds 
k£N Jo Jo Jri 

pT rs rsi p 

-16ImV/ / / / u(s2,x)Vu(s 2 ,x) 

^Jo Jo Jo JR d 

x V(4>ek)(x)dxdPk(s2)dsids. 

4. Blow-up for initial data with negative energy. The aim of this section 
is to generalize the well-known deterministic result stating that, for the 
critical or supercritical NLS equation, an initial data with negative energy 
yields a solution which forms a singularity in finite time. (See [20] , Chapter 
5, where more sophisticated results can also be found.) 



V{u(t)) = 



(3.5) 
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Theorem 4.1. Let uq, a and <f> satisfy the assumptions of Theorem 
2.1. Assume also that a>\, u £ L 2 (ft; E) n L 2<T+2 (0; L 2a+2 (R d )), f£ is a 
bounded function and, for some i > 0, 

(4.1) E(V(«o)) + 4E(G(«o))f + 8E( J ff(n ))t 2 + |t 3 m^E(M(u )) < 0, 



Remark 4.2. Clearly, if the energy of no is a.s. negative, then the con- 
clusion of Theorem 4.1 holds for some t > provided the noise is not too 
strong, that is, provided is small enough. 

Remark 4.3. By Sobolev's embedding theorem, the condition that f\ is 
bounded is fulfilled if is a 7-radonifying operator from L 2 (R d ;R) into 
W s ' p (M d ;]R) with (s - l)p > d. 

Remark 4.4. In [8], in the case of an additive noise, this result was 
proved under a severe assumption on a and the proof was rather technical. 
Under the assumptions made here, only a weaker result was proved. The 
argument used below can, in fact, be easily adapted to the additive case and 
the severe assumption on a can be removed. 

PROOF of Theorem 4.1. Assume that the conclusion of Theorem 4.1 
does not hold; then i < t*(uq) a.s. and we can take r = t as a stopping time 
in Corollary 3.5. In that case (3.5) can be simplified into 



where = 



then 



P(t*(u ) <t) >0. 



V{u{t)) = V(uq) + AG(u )i+ 8H{u )i 




(4.2) 






Let us set for t, r > 0, 



V{t, r) = V(u ) + AG(u )t + 8H(u )t 2 
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,2 -ad 



a + l Jo 
■t 



|^2ct+2 ds 



{t~ 8)111(8)11 

+ 4[\t-s) 2 [ \u(s,x)\ 2 fldxds 
Jo Jm d 

+ 4V / (t-s) [ \u(s,x)\ 2 x ■ V(4>e e )(x)dxd(3 e (s) 
^Jo JR d 

-8ImV / (t-s) 2 [ u(s,x)Vu(s,x) ■V(4>e i )(x)dxdl3i(s) 
eeN Jo jRd 

clearly, V(u(i )) = V(i, i ). Define the stopping time 

T k = inf{s E [0, t], \u(s)\ H i>k}, 
for k € N. Since the mass of u is constant, we have, assuming t < i , 

<"Tk 



E 



(i — s) / |it(s,x)| x ■ V(4>ee)(x) dx 



ds 



^m^El^J (t- syM(uo)V(u(s))ds 
^im^t 3 sup E(V(u(s AT k ))M(u )) 

sG[0,t] 

since = By Remark 3.4, we deduce that this quantity is finite. 

It follows that, when r = r^, the first stochastic integral above is square 
integrable and its expectation vanishes. Similarly, 



We 



Tk 



(t — s) 2 / u(s, x)Vu(s, x) ■ V(4>ei)(x) dx 



ds 



< 



lm^t 5 k 2 E(M(uo)) 



and, when r = r^, the expectation of the last term also vanishes. Moreover, 
using again the conservation of the mass, 

-t 



o 



(t-s) 2 / \u(s,x)\ 2 f£(x)dxds < p^m^M^o) 



and, recalling that a > |, we get that, for any fceN, 

E(V(t, Tfc)) < E(V(«o)) + 4E(G(«o))i + 8E(tf (u ))t 2 + |t 3 m^E(M(n )). 

We now choose t = i. By assumption, r^^i a.s. as /c goes to infinity. Thus, 
thanks to Fatou's lemma, we have 

E(V(u(t))) < E(V(uo)) + 4E(G(u ))t + 8E(H(u ))i 2 + |t 3 m^E(M(n )). 
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This last inequality contradicts (4.1), since V(u(i)) is nonnegative. Thus, 
the result is proved. □ 

The result stated in Theorem 4.1 is still valid if the expectation over f2 
in (4.1) is replaced by the expectation over an .Fo-measurable subset of £1. 
More precisely, defining, for M and H positive constants, 

(4.3) V&fi = {v£^, V(v) < M, G(v) < M, \v\\ 2 < M, H(v) < -&}, 

we can prove the following. 

Corollary 4.5. Let uq, a and 4> satisfy the assumptions of Theorem 
2.1. Assume also that a > ^, uq G S a.s. and f \ is a bounded function. Then 
for any M > and i>0, there is a constant H(t,M) > such that 

G Vm,h) > =^ Hr* (uq) < t) > 0. 

Proof. This is simply proved by taking H large enough so that 
M + AtM - 8i 2 H + ft 3 m M < 
and applying Theorem 4.1 with uq replaced by uq1q , where 

Note, indeed, that the solution of (2.2) satisfies u(t,uj) = for any t > if 
u(0,u)=0. □ 

5. Blow-up for any initial data in the supercritical case. In this section 
we assume that the nonlinearity is supercritical, a > ^. In that case we are 
able to strengthen considerably the result of Theorem 4.1. 

5.1. Main result. The main result of this article is the following. 

Theorem 5.1. Assume that | < a if d = 1,2 and ^<a<2ifd = 3. Let 
4> £ C^ 2 be such that ker^* = {0}. Then for any uq G S 2 with uq ^ and 
t > 0, we have 

P(r*(u )<t)>0, 

where t*(uo) is the existence time of the solution of (2.2) with initial data 
uq provided by Theorem 2.1. 

The proof of this result is given in Section 5.4. The idea is the same as 
in [8] where the case of an additive noise is treated. We choose t\, t2 > 
and consider a control problem where the noise is replaced by a control and 
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show that, for any initial data, there exists a control such that the solution 
U of the controlled nonlinear Schrodinger equation at time t\ satisfies 

(5.1) V(U(ti)) + 4G(U(ti))h + 8H(U(t!))4 + ^m^M (U (h)) < 0. 

Then, using the nondegeneracy of the noise, that is, ker^* = {0}, we know 
that the noise will be close to the control on the time interval [0,t\] with 
positive probability. Intuitively, this means that at time ti the solution of 

(2.2) will be close to U(t\) on a set of positive probability and that we will 
have P(u(£i) G g) > for some M and H as in Corollary 4.5 — the initial 
time being t\ instead of and i=t2- Thus, for the solution starting at 0, 
blow-up will occur with positive probability before the time t\ + £2- This is 
exactly the result since t\, £2 are arbitrary. 

In the case of a multiplicative noise considered here, this is much more 
difficult to justify rigorously. 

First, the control is a potential and it is not known whether the nonlinear 
Schrodinger equation is controlable by a potential. We prove in Section 5.2 
that, for any initial data, it is possible to construct a potential such that the 
solution at arbitrary time t\ has a very negative energy so that (5.1) holds. 
This construction only works in the supercritical case. It is not difficult to 
see that in the critical case, a = 4, such a result cannot be true. 

The second difficulty is that the Ito or Stratonovich products are not 
continuous with respect to the paths of the noise. Then, to replace the 
continuity argument used in [8] , we use more sophisticated probabilistic tools 
inspired by the characterization of the support of a diffusion (see [19]). We 
prove that, if the noise is nondegerate, the solution of the controlled problem 
is in the support of the diffusion associated with (2.2), which implies that 
with positive probability the solution will be close to the solution of the 
controlled problem. This is the aim of Section 5.3. Finally, we end the proof 
in Section 5.4. Technical results are proved in Section 6. 

The assumption on the smoothness of u$ is quite strong here since in the 
deterministic literature it is common to consider initial data in E. Our result 
can probably be extended to such data, however, this would considerably 
complicate the proof in Section 6.2 and, for clarity, we have decided to 
restrict our attention to this smaller class of initial data. 

5.2. Construction of a potential leading to blow-up. We first show that 
for a supercritical nonlinearity it is possible to build a control which acts 
as a potential in the nonlinear Schrodinger equation and has the property 
that, in a finite time, the solution has a very negative energy. Therefore, if 
the control is then switched off, and replaced by the noise, it follows from 
Corollary 4.5 that the solution will rapidly blow-up. 
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Proposition 5.2. Let u G S, u 7^ 0, Ti > 0. Let M = max(M(n ) + 
l;V(uo) + 4|G(uo)| + 1). Then, for any H > i/iere exists T2 < T\ and a 
potential f G L s (0,T 2 ;W 1 'P(R d )), for some s>l and some p with 1 < p < 
1 + swc/i t/iai t/ie solution of 

(5.2, ^-(AEZ + I^ + ZE/HO, 

t/(0) = u„ 

exists on [0,T2] and satisfies U{T2) G V a t/ h, where Vu R ^ s defined by (4.3). 
Moreover, if u G fT 2 (IR d ), i/ten / G C([0, T 2 ]; ff 2 ' 



Proof. We can assume that T\<\. Let -5 < a < a. We then consider 
the following supercritical nonlinear Schrodinger equation: 

i^-(AU + \\U\ 2 °U) = 0, 
(5.3) dt 

U(0)=u . 

The constant A > is taken such that the energy corresponding to equation 

(5.3) satisfies 

HaAu ) = f hvu (x)\ 2 - —^Mx)! 2 ^* 
jRrf Z Za + Z 

<-^(V(u ) + AT 1 \G(u )\). 

Since a < a, and a < 2 if d = 3, it follows from the deterministic theory 
(see [5]) that there exists a solution U in C([0, i*]; S) n L r (0, i*; W llP (M d )) 
for some t* > and (r,p) admissible with p = 2a + 2. Moreover, by the 
assumption above, 

^(«o) + mG(u ) + STfH^uo) < 

and, since § < a, this implies that J7 blows up at some time T* < T\. This 
is classical ([15, 20, 22]) and follows from the deterministic variance identity 
which states that if U exists on the time interval [0,Ti], then 

(5.4) V(U(t))<V(u ) + 4tG(u ) + 8t 2 H^ x (u y, t€[0,Ti]. 

But since V(U(t)) is nonnegative, this is clearly impossible. In fact, we have 

\U(t)\ H i — ► +00 and \U(t)\ L 2a+2 — > +00 wheni^T*. 
By the conservation of the L? norm and Holder's inequality, we have 

\U(t)\%t 2 +2 < 2A 2 ( 2 g + + 2 2) l^)l^ + c(X,a,a)\u \U 
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it follows 

H{U{t)) = 2^T2 mt)] ^ 2 " 2F+2 m) ^ + H **( U W 

~ ~W^ U{t) ^ + ^^c(X,a,a)\u \l 2 + H^(U(t)). 
Recalling that H^ t x(U(t)) does not depend on t, we deduce that 

H(U(t)) -oo when i^T*. 
Thus, there exists T 2 < T\ such that 

H{U{T 2 )) < -H. 

Moreover, by (5.4), 

V(U(T 2 )) < V(u ) + 4T 2 G(u ) + 8T 2 2 #*, A (^o) 

<y(n ) + 4|G(n )|<M 

since T 2 < 1 and H^ t x( u o) < 0. It is also classical that 

G(U(T 2 ))<G(u )+4T 2 H ajX (u ), 

so that G(U(T 2 )) < G{u ) < M. In addition, M{u{T 2 )) = M(u ) < M, hence, 
U(T 2 ) 6 V^h- It suffices to take f(t) = X\U(t)\ 2a - \U(t)\ 2a and U is then 
the solution of (5.2). 

Finally, we observe that if u € H 2 (R d ), it follows from [17] that U € 
C([0,T 2 ];H 2 (M d )). Now, if d = 3, f is not necessarily in G C([0, T 2 ]; if 2 (M d )), 
but it is easily seen that / G C([0,T 2 ]; W 1,p (K d )) for any p with 2<p< 
2d/(d — 2), and we may argue by density as follows: take /„ G C([0, T 2 ];H 2 (R d )) 
with f n ->■ / in C([0,T 2 ];W 1 'f(lR d )), where p < 2d/(d - 2) is such that for 
some g < 2d/ (d — 2), we have ^ = 1 — ~ = \ + \i that is, ^ + ^ = 2 • This 
is possible thanks to the fact that d < 4. Using then the method in [17], it 
is not difficult to see that the solution U n of (5.2), with / replaced by f n , 
exists on [0,T 2 ] and U n -> C/ in C([0,T 2 ];iI 1 (R d )). In addition, x£/ n xf7 
in C([0,T 2 ]; L 2 (M. d )), as follows from the same arguments as in [5], Corollary 
6.4.4. Hence, for n sufficiently large, f n fits the conclusion of Proposition 
5.2. □ 

5.3. The support result. We now show that, if the noise is nondegenerate, 
the solution U given by Proposition 5.2 is in the support of the diffusion 
associated to (2.2). More precisely, under some additional assumptions on 
the initial data, we prove that the support of the law of u(T 2 ) contains a 
neighborhood in E of U{T 2 ). In fact, following [19] and slightly generalizing 
our arguments, we could characterize exactly the support of the diffusion u. 
However, this is not the aim of this paper and, for clarity, we simply state 
and prove the result we need. 
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Proposition 5.3. Suppose that the assumptions of Theorem 2.1 hold, 
uq £ £ 2 , uq ^0, <p £ jC^' 2 and i/ie noise is nondegenerate: ker<^>* = {0}. Let 
T\ > 0, H > arbitrary, and T 2 < T\, U be given by Proposition 5.2; let u be 
the solution of (2.2) given by Theorem 2.1; then for arbitrary neighborhood 
V of U(T 2 ) in S, we foaue 

P(r*(u ) > T 2 and u(T 2 ) £ V) > 0. 

Proof. We adapt the method used to characterize the support of a 
diffusion and follow the same line as in [19]. 

We choose T > T 2 . Let (ek)k&n be a complete orthonormal system in 
L 2 (K d ;R) and (Pk)keN be a sequence of independent Brownian motions such 
that W = SfceN Pk<ft e k- We define W c = XifcgN ftkZk so that W c is a cylindrical 
Wiener process on L 2 (M d ;M). 

Since ker0* = {0}, we know that the range of (f> is dense in L 2 (IR c( ,]R) 
and thus in ff 2 (R d ,M). It follows that, for any n £ N, there exists g n £ 
C([0 1 T 2 ];L 2 (R d ,R)) such that 

( 5 - 5 ) 1/ _ < P9n\c([0,T 2 ];H 2 {R d )) - ~> 

where / is the potential given by Proposition 5.2. We set 

fn{t)=4>g n {t), t<T 2 , f n (t)=0, T 2 <t<T. 

We also assume that we have chosen (ek)keN m such a way that £ 
H 2 (R d ), k £ N. We denote by P n the orthogonal projector onto Sp(eo, . . . , e n ) 
and set At = -. We construct a piecewise constant approximation of the 
noise in the following way: 

<i(()= wm«!, * 6 [**,<*+ DM. 

and 

[we set W c {— At) = 0]. Then, we consider the following equation, for n £ N: 

(5.6) idu n - (An 11 + \u n \ 2u u n + / n n n ) dt = u n dW - dt - u n W n dt. 

By Girsanov's theorem (see [6], Theorem 10.14), 

W c ,n(t) = W c (t) - f\w c , n (s) - g n (s)) ds 
Jo 

is a cylindrical Wiener process for the probability measure 

dP n = D n d¥, 
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where 

D n — exp 

Jo 

Clearly, u n is the solution of 



7 \w C)n (s)-g n {s),dW c {s))-\ I \W c , n {s) - g n (s)\ 2 ds 



(5.7) idu 11 - {An 11 + \u n \ 2a u n )dt = u n dW n - l -u n F^dt, 

with W n = (j)W C}n . Since the law of the solution of the stochastic nonlinear 
Schrodinger equation does not depend on the probability space nor on the 
Wiener process, we deduce that the law of u n is the same as the law of u, 
the solution of (2.2). Therefore, for any Borelian set V, 

P(r*(u ) > T 2 and u{T 2 ) G V) = P„(r*(u ) > T 2 and u n (T 2 ) € V) 

= J ^-{t*(u )>T 2 and u n (T 2 )eV}D n dF, 

where we have denoted by t*(uq) the existence time for the solution of (5.7) 
with initial data uq. It follows 

P(r*(u ) > T 2 and u{T 2 ) € V) > 

if and only if P«(u ) > T 2 and u n {T 2 ) G V) > 0. 

We prove in Section 6.2 that 

(5.8) lim P(r*(no)<r 2 ) = 

n — >oo 

and 



n— >oo 



(5.9) Ve>0 lim P( t*(u ) >T 2 and sup \u n {t) - J7(t)| E > e ) = 0. 

t£[0,T 2 ] 



Therefore, for any V neighborhood in E of U(T 2 ), 

P(t* (u ) > T 2 and u n (T 2 ) GV)->1 when n -» oo, 
which clearly implies the result. □ 

5.4. Proof of Theorem 5.1. It is now easy to complete the argument 
sketched in Section 5.1. Indeed, for any t > and uq £ X 2 , we choose t±, 
t 2 > such that ii + 1 2 < t. We take M, as in Proposition 5.2, and choose 
H such that Corollary 4.5 is satisfied with M, and t = t 2 . 

Then we take T\ = t\ in Proposition 5.2 and obtain the control / and the 
solution U on the interval [0,T2] with T 2 <t\. It follows from the statement 
of Proposition 5.2 that U(T 2 ) € fj. 

We then use Proposition 5.3, with V = V ^ setting 

n T2 = {w G ^|t*(u ) > T 2 and u(T 2 ) € V}, 
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which is clearly an Tt 2 -measurable set, and which has positive probability 
according to Proposition 5.3, we deduce by Corollary 4.5, that the solution 
starting at time with initial data uiT^) will blow-up before t<i with a 
positive probability. But, if we denote this latter solution by u(t), then it 
is clear, by translating the Brownian motions, that the solution u of (2.2) 
satisifies for i > T2 : u(T2 + •) and u(-) have the same law. Hence, u blows up 
with positive probability before Ti + t2, and the result follows. 

6. Proofs of technical results. 

6.1. Proof of Proposition 3.2. The proof is based on a smoothing proce- 
dure, as in [10], and on the use of a truncated form of the variance: 

V £ (v)= [ e- £M2 \x\ 2 \v{x)\ 2 dx, veL 2 (R d ). 

The smoothing step is similar to the one in the proof of Proposition 3.1, 
which can be found in [10]. Thus, we omit it, keeping in mind that the 
computation below can be rigorously justified in this way. 

We apply Ito's formula to V £ {u{t)) (see [6], Theorem 4.17 for the form of 
Ito formula used here): 

dV e {u(t)) = {V^{u(t)), -iAu(t)) dt + (T/'(u(i)), -i\u{t)\ 2a u{t)) dt 

- ^(u(t)),u(t)F<t,)dt-(V^u(t)),iu(t)dW(t)) 

+ |Tr(y £ // Xu(t))(iu(t)<t>)(iu(t)(f>)*)dt. 

Note that, for v £ L 2 (M d ), 

(V^(u(t)),v) = 2Re(y d e- £lxl2 \x\ 2 u(t,x)v(x) dx^j . 

Hence, this term cancels if v is a real- valued function multiplied by iu(t). 
This is the case of the second and fourth terms above. 
We have 

(Ve(u(t)),-iAu(t))=4Jm I e- £lxl \l-e\x\ 2 )(x-Vu(t,x))u(t,x)dx, 

as follows after integrating by parts. Moreover, 
iTr(^'(u(t))(m(t)0)(m(t)0)*) 

= V / e- £lxl2 \x\ 2 (iu(t)<j)e k )(x)(-iu(t)(j)e k )(x)dx 
= / e" £|:c|2 |x| 2 |ti(t,x)| 2 F^(x)(i3; 
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and 

UV £ (u(t)),u(t)F <l> )dt = f e- £ W 2 \x\ 2 \u(t,x)\ 2 FJx)dx. 

It follows that the contributions of these two terms cancel. (This cancellation 
is natural. It reflects the fact that the noise is real and acts as a potential.) 
We obtain 

dV £ {u{t)) = Aim \ e~ £ ^ 2 {l- e\x\ 2 )(x-Vu(t,x))u(t,x)dxdt 

and for t < t*(uo), 

V £ (u(t)) = V £ (u ) 

+ 4Im / / e £|a:| (l-e|sc| ){x ■ Vu(s,x))u(s,x) dx ds. 
Jo Jm d 

In particular, if 

T k = M{te[0,T], \u(t)\ H x >k}, 
we have by the Cauchy-Schwarz inequality 

rtAr k 

V E (u(t A r fc )) < V £ (u ) + 4k / V}' 2 {u(s)) ds 



and by Gronwall's lemma 

V £ (u(t A T k )) < {4k 2 T + V £ (u ))e T . 
Letting e — > 0, this proves that, for t < r^, u(t) E £ and 
(6.2) V{u(t An)) < (Ak 2 T + V{u ))e T . 

Thus, u has trajectories in L oo (0,r/ c ; S) for any k G N. We now let k — > +oo 
and deduce that it has trajectories in L°°(0,r;S) for any stopping time 
r < r*(uo). It is now easy to let e — > in (6.1) and prove (3.4). 
Similarly, we apply Ito's formula to G(u(t)) and use the identity 

(G'(u(t)),-i(Au(t) + \u(t)\ 2r7 u(t))) 

= (2ix ■ Vu(t) + diu(t), -i(Au(t) + \u(t)\ 2a u(t))) 

= 2 f \Vu(t,x)\ 2 dx-^- f \u(t,x)\ 2c7+2 dx 

JR d <T + 1 JRd 

= AH(u(t)) + 2 -^\u(t)\%t 2 +2 . 
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We also have, after integration by parts, 
(G'(u(t)),-iudW(t)) 

= Rey^ / (2ix -Vu(t,x) + diu(t,x))iu(t,x)((f>ee)(x)dxd(3e(t) 

= V / \u(t,x)\ 2 x- V(4>et)(x)dxd0 e (t). 

Finally, a similar computation shows that 

(G'(u(t)), -iu(t)F^) + ^Tt(G" (u(t))(iu(t)4>)(iu(t)4>)*) = 0. 
Gathering these results and integrating in time yields (3.3). 

6.2. Proof of (5.8) and (5.9). We first introduce a truncated form of 
(5.6). We choose 3/2 < s < 2 and a cut-off function ■& G C§°(R) such that 

tf(z) = l, a;G [0,1], 

~d{x) G [0, 1], x>0, 

i?(a;) = 0, x>2. 

Then, for any i? > 1, we set 

We then consider the following equation: 

W«S - An^ eft - 6 R (u R )(\u R \ 2 °u R + eft 

= e R {u n R )u n R dw + Kr^ur) dt - e R {u R )u n R w n dt, 

with initial data 

u n R (0)=u . 

We have denoted by 
Note that, formally, 

e R (u n R )u n R dw + Ab,^(«s) ^ = e R {u R )u n R o div. 

It is not difficult to see that, for any R > 1, (6.3) has a unique solution with 
paths in C([0,T\; H 2 (R d )). 
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Step 1. For any R>1, there exists a constant C(R,T,cj)) such that, for 
any n € N, 

e( sup \u n R (t)\ 2 H2 )<C(R,T^). 

Vt€[0,T] / 

Proof. In order to lighten the notation in this step, we omit to explicit 
the dependence on R. All the constants appearing below are allowed to 
depend on R, T or </>, but not on n. We set 

u n = u% 9 = R , A = A Rj(j> . 

We use the Ito formula with F{u n ) = \u n \ 2 H2 = \ {I — A)u n j 2 and we have 



d\u n \ z H2 +2Re 



i / 9{u n ){{I- A) 2 u n )(\u n ru n + f n u n )dx 



dt 



= 2Re(-i J d 6(u n )((I - A) 2 u n )u n dW dx 
+ ZRe(-i J d ((I-A) 2 u n )A(u n )dx^jdt 

+ 2Re(i [e(u n )((I - A) 2 u n )u n W n dx) dt + 6 2 (u n ) £ \u n 4>e t f m dt 

We integrate between and t, and estimate each term of the resulting iden- 
tity. For the second term of the left-hand side, we write 



e( sup f 

\te[o,T]Jo 



<E 



T 



(u n )((I - A) 2 u n )(\u n \ 2a u n + f n u n ) dx ds 
6{u n ){I - A)u n {I - A){\u n \ 2a u n + f n u n ) dx 

,n\2a„.n i 



ds 



< cEU 9(u n )\u n \ H2 (\\u n \^u n \ H2 + \f n u n \ H 2)ds 

We have used an integration by parts and Cauchy-Schwarz's inequality. 
Now, we use the inequality 

|^li?2 < \v\l°° \w\h 2 + \ w\l°° \v\h 2 i v,w G H 2 (M d ), 

and the Sobolev embedding H s °(R d ) C L°°(R d ) to get that 

(6.4) I^H-ff 2 ^ \ v \h s o \w\h 2 + Mfl^o M.ff 2 > v, w S H 2 (M d ). 

This gives 

6{u n )\\u n \ 2u u n \ H 2 < cR 2u \u n \ H 2 < c\u n \ H 2 
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for another constant c which depends on R. By (5.5), f n is in C([0, T]; H 2 (M. d )) 
and its norm in this space is bounded independently of n. Therefore, by a 
similar argument as above, 

\fnU n \ H 2 <c\u n \ H 2. 

We deduce 



e( sup / 

\te[o,T] Jo 



(u n )((I - A) 2 u n ){\u n \ 2a u n + f n u n ) dx 



ds 



\u n \ 2 H 2 ds 



2 

In the same way, using (ft £ C 2 ' > we have 



and 



e( sup / 

\te[o,T] Jo 



E( sup / 

\te[o,T] Jo 



{{I-A) 2 u n )K{u n )dx 



£0V)K<^& 2 ds 

em 



\u n \ 2 H 2 ds 



\u n \ 2 H 2 ds 



Moreover, recalling that the noise is real valued and using a martingale 
inequality, 

Re/* / -i9(u n ){(I-A) 2 u n )u n dxdW 
Jo JR d 



E( sup 

ae[o,T] 



E sup 

ye [0,71 



El sup 

V*e[o,T] 



TlmV / d(u n )(I - A)u n (I - A){u n (/)e e )dxdPi 
Jo tm jRd 

x [u n A((j)e e ) + 2Vu n • V(0e/)] dxdfa 



2 \ V2 

x [n"A(^) + 2Vu n • V(0e<)] ) <ft j 



< C (E r^^d^i^Ki^i^i^ 



1/2 



\u n \ H 2\Vu n \ L3 \V(pe £ \ L6 y ds 
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by Holder's inequality. We now use the embeddings H s °(R d ) C L°°( 
ipo-i(Rd) c H i/2Q£d) c _^3( M d) and #i(R<i) c L 6 (R d ) to obtain 



E( sup 

.te[o,T] 



V)w n dxdW 



1/2 



Re A-0 [0(u n )((I-A) 

JO JR d 

< c( E A £ 0(u n )\u n \ 2 HH) \u n \ 2 H2 \<t>e t \ 2 H * dt 

<c+±E[ sup K|# 2 ). 

Vte[o,T] / 

It remains to estimate the most difficult term which we denote by A n (t). By 
similar manipulations as above, we have 

A n (t)= [ Re(i f 9{u n )((I - A) 2 u n )u n W n dx] ds 



Im( 

o V 



{u n )(I - A)u n {u n AW n + 2Vu n ■ VW n ) dx ds. 



We need to introduce a further notation. For s > and k £ N, we set Sk = 
kAt and Ug(s, s k -i) = S(s — s k -i)u n (sk-i)- Here (S(t)) te m is the group of 
operators associated to the linear Schrodinger equation, S(t) = e~ ltA , t G M. 
It is well known that, for any t, S(t) is an isometry on each Sobolev space 
H s (R d ), s£R. Now we split A n (t) as follows: 



-A n (t) 



Im / 9(u n (s))(I-A)(u n (s)-u n s (s,s k ^)) 

x (u n AW n + 2Vu n • VWn) da;) ds 

+ [\m([ 9(u n (s))(I-A)u%(s,s k ^) 
Jo \Jm d 

x (2V(u n (s) - ug(s, s fe _i)) • VW n ) dx) ds 
9(u n (s))(I-A)u n s (s,s k ^) 
x ((u n (s) - v%(s, s fe _ 1 ))AW n ) dx) 
(0K( S )) - 0(ug( S> Sfc -i)))(/ - A)«§(s, Sfc _0 
x (ug(s, s fc _i)AW n + 2Vug(s, s fc _i) • W n ) dx ) ds 



Im 



+ / Im 
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+ f Im( [ e(«S(s,s fc _i))(/-A)«S(a,a fc _i) 
Jo \Jw. d 

x (ug(s, Sk-JAWn + 2Vu^(s, s fe _i) • VVF n ) dx J ds, 



and in this expression, s^-i depends on s and is chosen to satisfy s € 
[ s k-i,$k[- We use similar arguments as above to deduce 

\A n (t)\<c f\\n n (s)-u n s (s,s k . 1 )\ m \W n \ H2 

+ \us(s,s k -i)\ H 2\u n (s) -Us(s,S k -x)\ H 'o\W n \ H 2 
+ \u n {s) -ug(s,S fe _i)|fl-« |«5(s,S fc _i)|H*o 

x |uS(a,afc_ 1 )| fl a|W r n |fl2)da 



o 



+ 



/' Im f / a fc _i))(7 - A)fi§(a, s k ^) 

JO \JB, d 



x («g(a, s fc _i)AW n + 2Vug(s, s fe _i) • VW n ) cfej 

= |Af(t)| + |^(t)|. 

The second term AV^ is a martingale and can be treated as the stochastic 
integral above. We obtain 



E sup A$(t) 
ae[o,T] 

< c(e [ T 2 {u n s {s, s k ^))\u n s (s, s fc _i)|^. s^OI^I^I^ dt] 



<c+±E sup |u n ' : 



4^ sup \u \ H 2 
te[0,T] 

The term A™ is more difficult to treat since EdW^I^) is of order (Ai)~ m / 2 
for any m G N. This factor is compensated by u n (s) — Ug(s, s k -i), which we 
now estimate. We write, using the integral form of (6.3) for u n , starting at 
* = Sfc-i, 

u n (s) -ug(s,s fc -i) 

5( S - T)[e{u n (r))\u n {T)\^u n {r) + /„(?») + A(/(r))] dr 

k-1 

+ if S( S -T)9(u n (T))u n (T)dW(T) 
Js k -i 
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S(s - T)9(u n (T))u n (T)W n dT. 

1 

Since H s °(M. d ) is an algebra, we obtain 
\u n (s) - ng(s,s fc _i)|^o 



< c(s - s k -i) + 



S{s-T)e{u n (T))u n (r)dW{r) 



Sfc-l 



+ c(s-s k . 1 )\W n ((k-l)At)\ H s 
and, for any m € N, 

(6.5) E(K( S ) - «g( a , Sfe _x)|^ ) < c(m)(s - s k ^) m ' 2 . 

Moreover, by (6.4), 

\u n {s) - U<j(s,S fe _i)|#2 

5(s-r)^(u n (r))n n (r)dVF(r) 



< c(s - Sfc_l)sup|li™|^2 + 
[0,s] 



+ c(a - Sfc_i) sup |u n |H2|Wn((fc - l)At)| H 2 
[o,«] 

and, by Cauchy-Schwarz's inequality, 



E(| U n ( S )-^( S , Sfc _ 1 )| H2 |W n (s)b2)<c(E(sup|« n |^ 2 



V[o,i 



1/2 



On the other hand, it is easily seen from (6.3) that, for any m G N, 



E sup |u"(t)lH»o ) <C(m,T,R)<c 



and we deduce 



E( sup Ai(t) I < c( El sup K(i)|^ 2 
-te[o,T] / V Vte[o,T] 



1/2 



It is now easy to gather all the preceding estimates and to conclude the 
proof of our claim. □ 

Step 2. For any R>1, there exists a constant C(R,T,ip) such that, for 
any n € N, 

e( sup \u n R (t)\l 2 ) <C(R,T,<f>). 
Vte[o,T] / 
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Proof. Thanks to the result of Step 1, it remains to estimate V 2 (u%(t)) = 
J Rd \x\ |tijj(s, x)| 2 dx. We use the ltd formula, recalling as in Section 6.1 that 
the computations can be justified by a smoothing argument. As above, we 
omit to write the dependence on R. We obtain 

dV 2 (u n ) = 21m / \x\' l u n Au n dxdt = -81m / \x\ 2 u n x ■ Vu n dxdt. 

We have, by Cauchy-Schwarz's inequality and integration by parts, 



\x\ 2 u n x ■ Vu n dx 



< {V 2 {u n )) 1/2 \x-Vu n \ 

< c(y 2 (u n )) 1/2 ((y 2 K)) 1/ V n l^ 2 + Kl 

<\V 2 (u n ) + c\u n \ 2 H2 . 



We deduce 



V 2 (u n (t))<e t 



V 2 (u ) + / \u 



,n\2 



ds 



We conclude thanks to Step 1. □ 



Step 3. The laws of {{u% W)) nm are tight in C([0, T]; XT') x C([0, T];H^ oc ( 
for any n < 2. 

Proof. We use the integral form of the equation satisfied by to 
get an estimate on its modulus of continuity. Again, we omit to write the 
dependence on R and all the constants appearing below do not depend on 
n. We have, for t\,t 2 £ [0,T] with t\ < t 2 , 

5(-t 1 ) u n (ti)-5(-i 2 )^(t 2 ) 

= -i [ 2 S{-s)[6{u n )\u n \ 2a u n + f n u n + A{u n )} ds 



S(-s)9(u n )u n dW(s) + i / S(-s)9(u n )u n W n ds. 



Using the fact that H s ° 
we deduce 



is an algebra and that S(—s) is an isometry, 



E(|5(-t 1 )^(t 1 )-S(-t 2 )u"(t 2 )|^ ; 

ft 2 



<c|t 2 -ti|+E 



ti 



S(-s)6(u n )u n W n ds 



H 3 o 



where again the constant c is allowed to depend on R. The last term is 
decomposed as A n in the first step into the sum of a martingale term 

S(-s k -i)9(v%(s, s fc -i))t$(s, s k -i)W n ds 



')) 
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and of a remaining term. It is not difficult to see that both have a second 
moment which can be majorized by a constant time \t 2 — ti\. Therefore, 

E(|S(-tiK(ti) - S(-t 2 )u n (t 2 )\ 2 H s ) < c\t 2 - h\. 

Now, since S(t) is an isometry and since, for v in H 2 (M. d ), we have 

(6.6) \(S(-h)-S(-t 2 ))v\ H s <c|ti-i 2 | 2 - S0 M H 2, 

as can be seen with the use of spatial Fourier transform and we easily obtain, 
thanks to Step 1, 

E(\u n (h) - U n (t 2 )\ 2 H s ) < C\t 2 - tl |2(2— 0). 

Thus, the sequence (u n ) neN is bounded in probability in C([0, T]; S 2 (M d )) n 
C 2 - s °([0,T];H S0 (R d )). This space is compactly embedded into C([0, T];T l r i(M d )) 
for rj < 2 and the tightness of (u n ) n ^ follows. 
We conclude thanks to the well-known identity 

E(\W(h) - W(t 2 )\ 2 H2 ) = \t x -t 2 \ H 2 ,2 
and the compactness of the embedding of H 2 (M. d ) into H^ oc (M. d ). □ 
Step 4. Passage to the limit. 

Let us fix a pair of subsequences (ip(n),^(n)) ne ^. We infer from Step 
3 and the Prokhorov and Skorokhod theorems that there is a subsequence 
of (u^ n \u^ n \w) which we still denote by the same letters, a probability 
space (Q,.F,P) and random variables (u™, u 2 , W n ), n £ N, (ui,u 2 ,W) with 
values in C( [0, T] ; S') x C( [0, T] ; S^) x C( [0, T] ; < c (l d )) for any rj < 2, such 
that for any n G N, 

£K,^,W") = £(4 W ,4 W ,W) 

and such that 

u^Uj as n->+oo, P a.s. in C([0,T];E") for r? < 2, j = l,2, 
PF n — > IF as ra->+oo, P a.s. in C([0,T];H^ QC (W i )) for r] < 2. 
Defining then 

.Ft = W(s), 0<s<t, i = 1,2} 

and 

J? = a{u](s),W n (s), 0< s <t, j = 1,2}, 

it is easily seen that IF and IF n are Wiener processes associated, respec- 
tively, with (Ff)t>o and (F™)t>o, with covariance operator (jxj>*. 
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It can be checked that for j = 1,2, we have 



(6.7) 



idu] - Au] dt - e R (u])(\u]\ 2a u] + f n v$) dt 



8 R (u])u] dW n + Ar^u]) dt - e R (u])u]W n dt 



and 



u](0)=u 



where W n is defined as W n but replacing W by W n . The difficulty now is 
to prove that the right-hand side of (6.7) goes to zero. In fact, we have the 
following result whose proof is very technical and is postponed to Step 6. 

Lemma 6.1. For j = 1,2 and any t G [0,T], 



in L 1 (£l,H s °(R d )). 

It is not difficult to let n — > oo in the left-hand side of (6.7) and, thanks 
to Lemma 6.1, we deduce that, for j = 1, 2, 



It can be seen that this equation has at most one solution so that the limits 
of u" and coincide: 



We deduce that for each pair of subsequences (ip(n),ijj(n)) ne ^, there is a 

subsequence of {u^ n \u R ^ n ^) which converges in law to a random variable 
(u\,u\). We now use the following elementary lemma, which was first used 
by Gyongy and Krylov in [16]: 

Lemma 6.2. Let Z n be a sequence of random elements in a Polish space 
E equipped with the Borel a-algebra. Then Z n converges in probability to 
an E-valued random element if and only if for every pair of subsequences 
(Ztp( n ) i Zif>(n)) i there is a subsequence of {Z v f n \,Z^t n \) which converges in 
law to a random element supported on the diagonal {(x,y) € E x E, x = y}. 

It follows that the sequence (u^) nS N converges in probability in C([0, T], T, 71 ) 
to a random variable u R . Moreover, proceeding as above, we can let n — > oo 
in (6.3) and see that u R satisfies the equation 

j idu R - Au R dt - R (u R )(\u R \ 2a u R + fu R ) dt = 0, 

1 ur(0) =u . 




idv,j — Av,j dt — 9 R (uj)(\uj\ 2a Uj + fuj) dt = 0, 
Uj(0) =u . 



U\ = U2- 
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Step 5. Conclusion. 

It is now important to emphasize the dependence on R. Let us define 
t r = min{t G [0,7], \u R (t)\ H s > R} and r£ = min{t G [0,T], |^(t)|^o > 
i?}. On [0, tr], ur satisfies the deterministic nonlinear Schrodinger equation 
(5.2) whose unique solution is given by U. Thus, if Ro = |t r |Loo(o,T 2 ;// s o(]R d )) 
and R> Ro, then tr > T 2 and ur = U on [0, T2]. Similarly, on [0, rjj], u R 
and u n , the solution of (5.6) coincide. Moreover, rjj and tr increase with R 
and 

Tr— >t* and Tr— >t*(uo) a.s. 

when — > 00, where i* is the existence time of U. It follows, since t* > 7%, 
that 

P«(u )<T 2 )^0 
as n — ► 00. Hence, (5.8) holds. Moreover, for any e < 1, 

P(t*(u ) >T 2 and sup \u n - U\v > e J 



<P (t*(u ) >7 2 and sup |u n -C/| s >e 
V te[o,T 2 AT% o+1 ] 

( sup |«JR 0+ l(t)-«B +l(t)| E > 

Vtefo.Tl 



< 



where we have used tr 0+ \ > T 2 . Now, this last term goes to as n goes to 
infinity by Step 4. We deduce that (5.9) holds. 

Step 6. Proof of Lemma 6.1. 

In order to lighten the notation, we omit writing the dependence on^i?. 
Also we write A instead of Ar^, u n instead of ffj and W n instead of W n . 
Thus, we have to prove 

I n {t)= f 6{u n )u n dW n {s) + f A{u n )ds- f 6(u n )u n W n ds^0, 
Jo Jo Jo 

when n — > 00 for any t G [0, 7] in L 1 (J7, H s °(M. d )). We split as the sum of 
L\ and J^, where 

J n(*) = E 0(u n (*)K(*) dW n {s) - f k+1 9{u n s {s, t k ^))u n s {s, t k -!)W n ds 
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and 



N-l 



+ 



tk+1 



J nW = E / WugMfc-OKCMfc-i) -e(u n ( s ))u n (s))w n ds 

k=0 ' 



A(u n (s))cZs. 



We have set i/% = kAt, k = 0, . . . , N — 1, tjy = T, where iV is the integer 
part of i/At, and t_i = 0. We have again used the notation itg(s,tfc_i) = 
5(s-t fc _i)ii n (t fc _i). Note that, for fc = 0, . . . ,N - 1, 



•fe+i 



where VK™ = <f> i W n is a cylindrical Wiener process. It follows 

N-l 



iUt) = E 

k=0 



tk ( ftk + l 



At 



(e(u n ( s ))u n (s) 

0(u s (T,t k .. 1 ))u s (T,t k -i))dr ) dW n (s) 



+ / e{u n {s))u n {s)dW n {s) 



N ~ l rt k 

+ E / 0{v n {s))u n {s)d( ( f ) {I-P n )Wl 

k=0 Jtk ~ 1 



Then, using similar arguments as in Step 1, 



N-l 

< E E 

k=0 



tk-l 



1 / ftk + 1 



At 



(e(u n (s))u n (s) 

- 0(us(t, tfc_i))itfif(r, tjfe-i)) c?t ) 



(6.8) 



+ / E(|^( U n ( S ))u n ( S )</»(/-P n )|>, S0 )d S 



In- 



E(\9(u n (s))u n (s)(l>\ 2 r o,s )ds 



<c(At)+c\cj)(I-P n )\ r0 ,s 
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and this latter term goes to as n — ► oo, since <p G C 2 ' S0 ■ 

Now, estimate (6.5) is not sufficient to handle i^(t). Let us define 

R(s,t k _ l ) = -i f S(s-T)[e(u n ){\u n \ 2 °u n + f n u n )-k(u n )]dT 
Jtk-1 

-if S(s-T){8(u n )u n dW n (T)} 

Jtk-1 



+ i / e(u n s (s,t k . l ))u n s (s,t k . 1 )dW n (r) 
+ i f S{s-T)[6{u n )u n W n dT] 

-if 0(uS(s ) t fc _l))ti§(s,t fc _l)W' n dT 



so that 

u n (*)-T#( S> tfc-i) 



2?(Mfc-i)-*0(«!KMfc-i))u.KM*-i) f <W n (r) 

Jtk-i 



+ ie(u n s (s,t k _ 1 ))u n s (s,t k _ 1 ) / W„dr. 

We can prove, using Cauchy-Schwarz's inequality, Step 1, (6.5) and (6.6), 
that for any 5 < 2, 

E(\R{s,t k - 1 )\U)<cA^' 0+1/2) . 

It follows 

(6.9) E(|i?( S ,t fe _ 1 )|^o|Wn|^o)<cAi 2 - so . 
Next, we rewrite as 



k=0 



J („.n i 



a' g (u n s (s,t k ^)) ■ (u n s (s,t k ^)-u n (s)))W n d S 
N-l r t k+1 

£ / Jgivfifatk-!)) ■ R(s,t k -i)W n ds 

7~ r> JU 



k=0 
N-l 



+ ]T [ (Mu n (s)) - A(u n s (s,t k ^)))ds 
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N-l r t k+1 

Y / 0(u n s (s,t k ^))a' e (u n s (s,t k ^)) 
f~1, Ju 

'tt$(M*-i) f dW n (T))w n ds 



k=0 Jtk 



N-l . t 



k=0 



mS(a,t fc _i) jT W n dTjW n ds 

+ Y / A(«g(s,t fc _i))da, 
fc=0 Jtk 

where cr' s (v) ■ w is the differential in the direction w of the mapping v i— » 6(v)v 
which is C°° from H s ° into itself with all derivatives bounded. By the Taylor 
formula and (6.5), the first term clearly goes to zero in L 1 (fi; H s °) as (At) 1 / 2 . 
This is also the case for the third term. For the second term, we use (6.9) 
and see that it goes to zero in L^fijfF ) as (At) 2 ~ s °. 

Let us denote by J(t) the sum of the remaining terms. Using that 

A(u) = \0(u) Y a 'e( u ) ■ (-iu4>ee)4>ee, 

we get 

j(t) = Y / tfc-i)) 

1- n Jtk 



k=0 



^(uS(s,tfc-i)) 

iu]K«,*fc-i)(T dW n (r)- f W n d T ))w n 
\Jtk-i Jt k -\ J J 

+ lY a 9( u s( 3 > t k-i)) • (-ius(s,t k -i)<i>ei)(f)et 



ds. 



Now, recalling the definition of 9, we have 
and we can write 

J(t) = J 1 (t) + J 2 (t), 
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with 



Ar ~ 1 rt k +i 

Ji(t)=J2 / i0 2 («3( s ,t fc _i))u§(*,t fc _i) 



k=0 



dW n (r) 



W n dT)W n -±F 4 



ds 



and 



o N ~ 1 rt k+1 

J 2W = ^E / 0(«§(M*-i))*' 

« fc=0 J<fc 



,/|ug(s,t fc _i] 



Lff s o 



i? 2 



ttg(s,tfc_l),iUs(s,tfc_i) 



tfc-l / / fPo 



X ^(s,tjfc_i)ds. 



Writing 



*fc-i 



dW n (r) 



tfe-i / 



tfe-i 



tfc-i 



it is then the result of tedious — but standard — computations, based on the 
independence of the increments of the Wiener process, that each bracket 
in the right-hand side above gives, when inserted into Ji(t), a term which 
goes to in L 2 (0; L 2 (R d )), hence, also in L 2 (0; H s °(R d ))— since it is clearly 
bounded in L 2 (f}; H 2 (R d )). 

It can be proved in the same way that J2W goes to in L 2 (Q; H S0 (M. )) 
as At goes to (or n goes to infinity). 

Therefore, I 2 (t) goes to zero in L 1 (f];i7 s °) and, recalling (6.8), Lemma 
6.1 is proved. 
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